THE CAMBRIDGE 


MATHEMATICAL J OURNAL. 


IL—ON THE PROPAGATION OF A WAVE IN AN 
ELASTIC MEDIUM. 


THE following investigation of the equation for the propagation of 
a wave in an elastic medium, may be considered as supplemental to 
the two articles on the Undulatory Theory in the first and second 
Numbers of this Journal: though perhaps it ought more properly 
to have preceded them. But as we do not pretend to offer a com- 
plete treatise on the subject, the order of the articles i is not a matter | 
of much consequence. 


1. Light, according to the wave theory, consists of undulations 
propagated in an elastic medium. From the laws of the interference 
of polarized light, it appears that the constituent vibrations are 
transverse to the direction of propagation ; that polarized light con- 
sists of vibrations in one direction, or perpendicular to one plane, 
which in Fresnel’s theory is the plane of polarization, and in 
Cauchy’s theory is perpendicular to the plane of polarization ; and 
thirdly, that there are no vibrations, or at least none capable of 
producing the sensation of light, in the direction of propagation. 
Fresnel assigns as the reason for this, that the resistance to com-— 
pression of the medium is so great, that if a disturbance takes place 
motion is propagated to a great distance before the disturbed mole- 
cule returns to rest. Hence, a wave of great breadth will be pro- 
pagated, but of feeble intensity, as the vis viva remains constant. 

2. To this we may add, that the possibility of the propagation 
of motion to a distance depends upon this, that waves of constant 
breadth can be propagated: if we examine our equations, we see 
that this depends upon the disturbance being a function of vt—2z 


. 
; 
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and the partial differential equation which gives this form of solu- 
tion is obtained, in the simplest case, namely that of motion in one 
direction, by supposing ‘the elastic force proportional to the con- 
densation, and the condensations to be infinitesimal. Without 
these restrictions the equations cannot be solved; but from several 
considerations it seems probable, that a particle once disturbed does 
not immediately come to rest, and that the waves increase in 
breadth. We may see examples of this in the waves propagated 
when a stone is thrown into still water. -But in water the velocity 
of a wave depends on its breadth; and thus, though the breadths of 
_ the waves increase, they do not interfere with each other. 


3. It appears probable, from the change in the character of 


sounds caused by distance, the loss of sharpness and distinctness, ae 


that the breadth of a wave of air increases. If we suppose the 
same to be the case with those waves of light which consist of vi- 
brations in the direction of propagation, but to a greater degree 
than in sound, we shall have, in the first place, a rapid diminution 
in the vis viva of. each wave ; and also, as their velocities are nearly 
independent of their breadths, each wave will overtake its prede- 
cessor and be overtaken by its successor, and there will be an 


interference which = speedily obliterate all traces of the disturb- 
ance. | 


4. We have now to consider that complex system of disturbances 
which constitutes an undulation. 


We shall suppose the molecules of the media which we consider 
to be so arranged, that the three rectangular directions, or the three 
axes of elasticity, have the same direction at every point. These 
we shall take for coordinate axes, and shall suppose that the elas- 
ticity called into play is greatest for a displacement parallel to the 

axis of a, and least for one parallel to the axis of z, or that the axes 
of x, y, z are the axes of greatest, mean, and least elasticity. 


5. If we suppose a range of particles, in a plane parallel to that 
of yz, to be moved upwards, or parallel to the axis of z, the 
resultant of the attraction of all the others will tend to pull it back 
along the same direction with a force proportional to the displace- 
ment: the same thing will happen if it receive a displacement pa- 
rallel to the axis of y, but for the same displacement the force of 
restitution will be greater than before, as we have supposed the 
elasticity for displacements parallel to y to be greater than for those _ 
parallel to z. Let a?, b%, c? be the coefficients of elasticity, so that — 
in the first case the force tending to bring back the range of parti- 
cles is c? X displacement, and in the other 62% displacement. If the 
influence of one range of particles does not extend beyond the 
nearest range, the force of each of the iy ape ranges will pull the 


displaced range back with a force = 5 X displacement. 


Let us now consider a succession of ranges to be displaced by 
different quantities, all parallel to the axis of z, the rauges being all 
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at equal intervals. Let the range at a distance 2 from the ortirin 
be displaced by a quantity z,, that at a distance a +A by a quan- 
tity z,,,, and so on. The range at x will be pulled up by the 
range ata -+h with a force (z,,, — z,) Az,, and down 
by the range at a —h witha “force 4 c°Az,_,: thus from these two 
it will be pulled wp with a force equal to their difference, or 
2 £0? (Az, - Az |: | 

and if M be the mass of the range at a, we have for its motion 

6. The general solution of this equation of mixed partial differ- 
ences has been given in Article II. of our last Number; and there it 
is shown, that the velocity of propagation depends on the length of 
the wave, and thus the phenomenon of dispersion is accounted for. | 


In the present case we shall take the approximate solution, where h 
2 


| dz 
is very small. By for Az 


, our equation becomes 


rah 


This is the well-known equation, which shows that a wave will 


h 
be propagated along the axis of x with a velocity ¢ ET If the 


vibrations take place parallel to the axis of y, there will be a wave 


propagated with a velocity 6 Fresnel states, that he has 


h 
Vv 2M 
satisfied himself, by experiment, that the velocity of a plane wave 
depends solely on the direction of the constituent vibrations, and 
not on the direction of the wave. It would appear from this, that 
the factors of 6 and ¢ in the above expressions must be constant for 
different positions of the wave; and it results, that plane waves, 
consisting of vibrations parallel to the axes of 2, y, or z, will be 
propagated perpendicularly to themselves with velocities which are 
-proportional to a, b, c, or to the square roots of the coefficients 
of the elastic force for the disturbance of a single range. 


7: If the range of particles in the plane of yz receive a dis- 
turbance in that plane, but not parallel to one of the axes, we must 
resolve this displacement into two, one parallel to each axis, and 
there will result ¢wo waves, propagated with velocities 6 and c re- 
spectively along the axis of x, and corresponding waves propagated 
in a negative direction. 


And this result follows, whether the displacement takes place 
along a line, as in polarized light, or in some more complicated 
manner, as in common light; and thus any kind of disturbance in 
one of the coordinate planes will produce two waves, polarized at 
right angles to each other, and of different refrangibility. 
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8. We have now to consider the general case of. a plane wave 
inclined at any angles to the coordinate planes. A vibration in this 
plane will produce a force of restitution, which will have a different | 
direction from the displacement. Resolve this force into two, one 
along the direction of displacement, the other perpendicular to it; 
if this last portion be also perpendicular to the front of the wave, it © 
- may be neglected, and a wave will be propagated with a velocity 

proportional to the square root of the elasticity resolved along the 
direction of vibration ; if it be not perpendicular to the front of the 
wave, it cannot be neglected. We have, therefore, first to inquire 
whether this is the case for any direction of vibration in the plane ; 
and this leads us into the subject which was discussed in our last 
Number. 


A. S: 


I]. —CIRCULAR SECTIONS IN SURF ACES OF THE 
SECOND ORDER. 


In determining the circular sections in the Surface of Elasticity, 
Fresnel has made use of a method which is very readily applicable 
to surfaces of the second order; and as it has not yet been intro-_ 
duced into any work on Analytical Geometry, it may be useful to — 
insert it here. 


Taking first the surfaces which have a centre, let their equation 
be 


and let this be cut by a atid 


which we suppose to pass ial the centre, as all sections made 
by parallel planes are similar. Let this plane also cut the sphere 


Now, as m, ”, r are indeterminate, we can so assume the position 
of the plane and the magnitude of the sphere, that the circular sec- 
tion of the surface (1), if it exist, shall coincide with the section of 
the sphere ; and if these coincide, the equations to their projections 
on the plane of zy must be identical, which gives us conditions for 
determining m and 2. Substituting for z in (1) and (3) its value 
from (2), we get | 

(P + 224. (P’ 4+ P’n?) y? + 2P"mnay H......(4), 

(1 4+ + (1 4+ 4+ 2mnazy =7"......(5). 
Comparing each term separately, those involving xy will coincide 
if either | 

H 


? 
; 
i 
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Taking the first condition and comparing _ other terms, 


| r, and H pee 
which gives + P"n? = P(1 4 7n?), 


and a= + Y 


If we suppose n = 0, we find in the same manner 
P — P 


The third condition leads to no result, and therefore i is not to be 
considered. 


In the ellipsoid, Yr, wea positive, and 


This shows that the value of x is impossible, and that of m possible ; 
therefore there are two directions arising from the doubtful sign 


in which the ellipsoid may be cut in cireular sections, determined 
by the equation to the cutting plane, 


pe 


In the hyperboloid of one sheet P” is negative, and the value 
of m is possible and m impossible. In the hyperboloid of two 
sheets P’ and P” are both negative, and. P” < P’, m is possible _ 
and 2 impossible. It is true, that for a plane passing through the 
centre the section is impossible, but a plane drawn parallel to this 


at a sufficient distance from the centre will cut the surface in a 
circle. 


m= 


The equation to the surfaces without a centre is 


py? + = ppr. 
Let this be cut by a plane 
mz NY, 
which also cuts the sphere | 
ayy? 4+ 227 | 
The equations to the projections of the sections on the plane of 
zy are | 
py + p2— mpp'z — npply = 0, 
(1 n?) y? + (1 + m*) 2? 4+ Qmnzy — — QZnry = 0. 
In order that these may coincide, the term involving zy must 
vanish, which will be the case if m =0 or n=0. 


If m = 0, then 1 +n? andn=+ 
Pp 


If n = 0, then 1 and m= + \/ 
P 


‘ 
| 
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In the elliptic paraboloid p and p’ are both positive, and accord- 
ing as p' is greater or less than p, the first or second is to be taken, © 
the other becoming impossible. In either case there are two series 
of circular sections corresponding to the positive and negative sign. 
In the hy perbolic paraboloid p or p is negative, so that there are 
no sections in which it is cut ina circle. This would appear also 
from the nature of the surface, as it can never be cut by a pane in 
a. closed curve. 
The same method | may be sipplied to the oblique cone, so as to- 
_ determine the sub-contrary sections. 
By determining the circular sections in this manner, it is seen at 
once that any two belonging to different series are situate on the 
same sphere. 


F. 


II.—PROPOSITIONS IN THE THEORY OF NUMBERS. 


_ Tue proposition, “The continued product of any m consecutive — 
integers is divisible by 1.2.3...m,” has been strictly proved by 
induction ; but this mode of proof only shows, if we may so speak, 
that the proposition must be true, without showing why it is true. 

We have seen it reasoned, that since the product of any m con- 
secutive integers is divisible by each of the natural numbers up to 
m separately, it is divisible by their product: but a little consider- 
ation will show that this does not follow. | 

The problem from which the above proposition is here deduced, 
has been proposed and solved before, but not, that we are aware, 
in the same form. 


Prop. 1. If p be a prime’ oii the index of the highest 
t phy of p that is a divisor of the continued product 1.2.3...M, is 

— §(M) 

p-—i 
expressed in the scale, whose radix is p. 


We shall adopt the notation used by several Continental mathe- 
maticians, M! for 1.2.3...M. 


Let the quotient of M by p be M,, and the remainder a,, then 
M, P+ % | 
and the multiples of p in the series 1, 2, 3, &c. M, are 
p, 2p, 3p, &c. Mip; 


wherefore, if K denote the product of the remaaning numbers in 
the former series, 


p-2p.3p...Mip, 
M! 


; where S(M) represents the sum of the digits of M 


| 

‘ 
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Let r be the index of the greatest gene of p in Mi, 
then from the | 
| r= M, + 
For the same reason, .if 

M,= M,p+a4,, and M,=M,p + a,, &c. 
and 7, Ts; ie. be corresponding quantities for M,!, M,!, &c. 


My Ty 


M,, being that quotient less than p at which we must t finally arrive 
Adding all these ew 


= M, M, eee eee $+ M 
But, adding the equations | | 
= Mp + 


= + a, 
M,, a3 


and observing that a,, a,, &c. a, are the successive digits of M ex- 
pressed in the scale of p, we have 3 


M + M, + M,+...+ M;,= (M,+ M,+.. + M,) p+ S(M): 
whence M, + M,+ . M,, 
and the proposition is proved. 


The following proposition is also useful. 


Prop. 2. If M and N be two numbers, of which M is the 
ora and S(M), S(N) the sums of their digits to any radix " 
S (M — N) is not less than S(M) — 8 (N).. 


For let M = '+....4 4,7" + 
then | 
M— N=ay? +4, —b, + (a,_,— +... 
Now, as long as the digits of N are not greater than the corres- 
ponding digits of M, the digits of M — N will be the excesses of 
those of M over those of N. But suppose 6, greater than a,, then 


we form the digits of M—N by writing that part of the difference 
thus, 


e 
44 
3 
* 


2 
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and the sum of these two digits is 
hence S(M — N) may be greater, but cannot be less, than 
S(M) — S(N). | 

Prop. 3. The product of any m consecutive integers is divisible 
by 1.2.3... m. | 

By Prop. 1, the index of p in Nii and in (N-+m)!, 
+ m — S(N + ; hence the index of p in 

(N+ 1) (N42). (8 +m) or 
is the difference of these quantities, or 
m —S(N 4 m) + S(N) 
| 

m—S§ (m) 


The index of p in m! is 3 and, by Prop. 2, ‘Making 


M = N + m, S(m) is not < S(N + m) — S(N), 
therefore m — S(m) is not > m—S(N 4+ m) + S(N);_ 


and the index of p in m! is not greater than the index of Pp in 


(N +1)(N 4+ 2)...(N +m). Hence, if 
(N 4 1) 4 2)..0N + m) = 
ais not > a, > b, y not &e. 
(N 4+ 1)(N 2).. is divisible Py 203. 


8. S. G. 


IV, —NOTES ‘ON FOURIER’S HEAT. 


. THe method employed by ears to integrate the partial differen- 
tial equations which occur in the Theory of Heat, is to assume 
some simple form of a singular solution, and afterwards to extend it 


go as to include all the” circumstances of the problem. It is in ef- 


fecting this that he has displayed the great resources of his analysis, — 
and imparted so great an interest to his work by the variety and 
ingenuity of his methods. Indeed there is a freshness and origi- — 
nality in the writings of Fourier which make them in no ordinary 
degree arrest the attention of the reader. But however much we 
may admire the means by which Fourier has overcome the diffi- 
culties of the problems he had to deal with, yet it seems more 


> 
mre 
be 
WN. 
‘ 
‘ 
> 
é 
a 
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agreeable to the usual style of mathematical investigation to deduce 
a result by limiting the general solution by means of the conditions 
of the problem, than by extending a particular case. 


That this may be sometimes done with even more readiness than 
by Fourier’s method, will be seen by the following solution of a 
problem given in p. 161 of the Théorie de la Chaleur. We may 
_ remark, that there is in general no difficulty in the solution of the 
partial differential equations, but only in the proper determination. 
of the arbitrary functions in the —— so as to suit the con- 
— ditions of the problem. 
If a rectangular plate, bounded by two infinite parallel edges, : 
have one of its extremities kept at a constant temperature 1, while © 
- the infinite edges perpendicular to the heated edge are retained at 
a constant temperature 0, the equation from which the temperature 
is to be determined is 


dx? 4+ "7,2 = () (1), 


where v is the temperature at the point Xs Y; the origin being at the | 
middle ‘point of the heated edge, the axis of a bisecting the plate, 
and the axis of y parallel to the heated edge. For the sake of 
shortness Fourier represents the breadth of the plate by z. : 


~The solution of the equation (1) by the method of the separation 
of the symbols of operation from those of quantity, is 


v = COS (y 5 (x) + sin (y (2) (2), 


(2) and W (ar) being functions of x. And it may also 
be put under the form | 

where F(a) = + and = (2) — 


Now, on Ledlind at the circumstances of the problem, it will be | 
seen that it must be subject to the following conditions : 


1st, » must be symmetrical with regard to y and — y. 


gna, y = 0 when y = 5 or — oo whatever x may be. 


3rd, » = 1] when 2 = 0, whatever 7 y may be. 

4th, » must be very small when z is very large. 

From the first condition we must have (a) = 0, as otherwise 
the second term would change its sign when —y is put for y. 
Hence we have only 


=) (2) 
(3) 
By the second condition, putting 5 for y in equation (3), we have 
ad | 
O = cos{- — ces 
(FF) (2) (4) 
. 
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Now this is in fact a linear differential equation with constant 
coefficients, and of an infinite order. By the principles laid down 


in Art. V. of our first Number, we can integrate this equation if we 


know the roots of the equation cos ¢ z) = 0. Now these are 


+1, +3, +5, &e. 
being in number infinite. Hence the solution of (4) is 
? (x) = 32 5x 
+ 4+ 4 &e. 
the number of ‘aes and tei constants being infinite. By 
the fourth condition it appears that the second line of (5) must 


disappear, as otherwise v would be very i when 2 is very large. 
Hence we must have 


C',=0, C,=0, C,=0, &e. 
and equation (5) is reduced to | | | 


~~ and vw becomes 


By the third condition v =1 whenz= 0. If then we expand the 
symbol of operation in (7), operate on each term separately, and 
then make x = 0, we shall find 


1 = C, cos y + C, cos 3y + C, cos 5y + &e......(8); 


g ond 


In order to determine the arbitrary constants, we shall follow | 
Fourier's method of definite integrals. If we multiply both sides 


where y is contained between the limits — 


of (8) by cosy dy, and integrate between the limits +5 and — 5 ; 


all the terms except the first will disappear, as they can each be 
decomposed into the cosines of even multiples of y, which, on 
Integration, vanish at both limits. Hence we have » 


Coe ¥ = cos? y = (1 + cos 2y), 
whence we find C, = i. 


In a similar manner we should find 


4 


4: 
3 Cs 


— 


and so on. 


. 
> 
: 
Ad 
‘ 
4 
« 
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| Substituting these values in equation (7), it becomes 


| 
Now if we expand the sign of operation, and apply it to sucha | 


term as e~"*, we shall find that it becomes cos nye~"*. Hence the 
expression for v becomes 
T 


which is one form of the solution whieh Fourier gives. ‘Tt may 
_ easily be reduced to a more simple form. For if we substitute for 


the cosines their we have 


2 (ety =I) 4 4 _ &e. 
which, by Gregorie’s series, become | 


| 2cosy \ 
which i is the simplest form that the ‘expression can assume. 


V._DEMONSTRATIONS OF SOME PROPERTIES OF _ 
SURFACES OF THE SECOND ORDER. 


To the Editor of the Cambridge Mathematical Journal. 


Srr,—I observe the author of a Paper in your first Number 
says, he is not aware that any person has made use of the 
symmetrical form of the equations to astraight line: he will find 
them employed in Cauchy’s “ Lecgons sur les Applications du Calcul 
Infinitésimal.” 

Perhaps the following solutions of two common problems may be 
acceptable to some of your readers. 

1. To find the locus of the bisections of — chords in the 
surface whose equation is 

Az? + By? + C2? — 


v = tans + | 
— — 
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_ Suppose any one of the chords meets the surface in the points — 
_ whose coordinates are 2 Yy 2» Tp Yy» 2,3 then if a, 3, y are the 
angles which the chord makes with the three axes, we have 


— 


Let a, y, z be the coordinates of the middle point of the chord ; then. 


therefore, combining these equations with (1), we have 
y2—y? 22-22 
xreosa  ycosp  Zzcosy 
Now, from the equation to the surface we have | 
A —. 2?) + B (Yo — ¥?) + — = 0, | | 
which, combined with (2), gives immediately | 


Now suppose & », Z are three semiconjugate diameters of the 
surface proposed, we have of course _ 


+ By y, + Cz,2z, = 0, 
Az,t, + Byy, + Cz,z, = 0, 
Az,2, + By,yi + Cz,z, = 0; 


Ax cosa + By cos 6 4 Cz cos y = 0, 
which is the equation to the locus required. : 
' 9. To find the locus of the intersection of tangent planes drawn : 
parallel to any system of conjugate diameters in the same surface. 
I shall first observe, that the projection of the diagonal of a 
parallelogram on any line whatever, is evidently equal to the sum 
of the projections of any two contiguous sides on the same line. 
Whence it easily follows, that the projection of the diagonal of a 
: parallelopiped is equal to the sum of the projections of any three of 
its edges, of which no two are parallel to each other. | | 
Hence, if , n, are three straight lines drawn from the origin to. i 
three points P,, P, Pa; whose coordinates are 
Loy Yo» Z, and if through the points P,, P,, P, we draw planes a 
parallel respectively to the planes of nZ, &, &, these three planes 
will intersect in some point P, and a line drawn from the origin to 
P will be the diagonal of a parallelopiped, of which é, 7, Z are three 
edges. And if we call a, y, z the coordinates of P, and project the 
diagonal and the three lines &, n, ¢ on the three axes successively, we 
have evidently, by the principle above mentioned, | 
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and if we multiply these equations by 2, and add them to the 
three following, viz. 


An? + By? + Cz? = 

Aa,’ + By? =f, 

An? + By, + Cz? = 
we obtain immediately 


A (a) +2,+2 + C(z, = 


or r by (3) 
Aa? + By? + Cz? — 3h, 
which is evidently the equation to the locus required. 


I do not know whether this solution is. new, but I do not re- 
collect to have seen it before. 


The equations (3) lead immediately to the common formule for 
the transformation of coordinates; for if we call a, By yo, ay By yp 
Go, Io) Yo» the angles which ¢, 7, ¢ make with the three axes respec- © 
tively, we have x)= cosa, y,= — By =F Cos yy, &c. and 
the equations (3) become. 

x= cosa, + n COS a, + ¢ COs a,, 

y = + 7 cos /3, + cos 

E COS 4. n COS y, + COS 
where é, n, ¢ are evidently the coordinates of P, referred to lines 
drawn from the origin through P,, P), P, as axes. 


I am, Sir, &e. M.N.N. 
Oxford, 31, 1838. | 


VI.—ON GENERAL DIFFERENTIATION.—No. IL. 


In our first Number we investigated the general differential coefhi- 
cients of e"* and x”. We shall now proceed to find those of other 
simple functions, and to give some examples of the: application of | 
the 

- We found that our formule gave an infinite value for — 


3 


when 2 is not, and m—a is, a positive integer. But by using the | 
complementary function, a finite value of another form may be 


obtained, in the same manner as the value of fz — may be derived 


from that of fx"dz. 


Let a —n in formula (N), and » — a in formula (O), (p. 19), be 
assumed equal to t+, when 7 is an integer, then the second 


a 
‘ 
& 


4 
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. +2 
member of each will be a multiple of ni » which is 
sin (7 + 


‘infinite when B =0, but by sagas in a part of the complementary 


x" 
sin 
log when B=0. “Hence i in these cases (N) and 


function becomes 
(-1)" 


( 0) become respectively 
i 


which i is fraction whose | 


value is - 


d* log 


2. The value of 


may be readily derived from the 


formula for the differential coefficients of powers of 3 x. _ For since 
Hence, if a be positive, : 
d* | 


us 


| 

3. If it be known that | 
= ¥(2), 


it is evident that 


hence the general differential coefficients of all rational functions of 


_ a can be found, by preparing them as for integration. Also, if the 


differential coefficient of a function to any integral index be a 
rational function, its general differential coefficient may be found. 
As an example, let 

y = 2, 


therefore, using the formula (P),, | 
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Substituting for a its value tan y, observing that 
v—1 = cos (2r + 4) m7 + /— 1 sin (2r + 4)z, 
and reducing by De Moivre’s theorem, we obtain | 


a—l 
); is that of x, is 


P(- 
The value of 


when a is positive, and 


| ( ) sin (— an). (1— a) 
when a is negative. The formula fails when a is a negative integer. 


4. Let us ascertain in what cases the differential coefficient of a 
constant — is not zero. 


Since C = C2", | 


which will not be zero if P(- a) be not infinite, which is only 
when — a is a positive integer, or a a negative aaa, that is, in 
_ the case of common integration. 


5. The differential coefficient of infinity may be finite. For 
«© = CP(m) x", if m be not a positive integer: therefore 


P(m) P(n) ‘ 
which will be finite if m be a positive integer, and x — a not. 


6. We proceed to find the general differential coefficients of 
cos and sin x. 


Since cos x + 1 ta 


= feos (2r +1)7 4 sin (2r + 4) +V— 1 sin 2) 


In like manner, 


d* cos x d* x 


= cos {(2r' 4+ + sin + 4) an + 2}, 


r’ being some integer, not necessarily the same as r. Adding and 
subtracting, | 


| 

3 

* 
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d* cos & 
= feos — + sin — 
COS 3(7 + 7 ar v 
sin $(7 + 7) + 4) ar + 2, J 


Since r and 1 r are both arbitrary, there is no relation between 
r +r’ and r —7, except that they must be both odd or- both even. 
If abe a rational fraction whose denominator, when in its lowest 
terms, is 2, r+ 7’ must go through all the 2n values, — 7 


before the variety of values of the trigonometrical expression into — 
which it enters can be exhausted; and the same is true of r — 7’. 
Now any of the m odd values of 7 + 7° may be combined with 
any of the v. odd values of r—7’, and thus x® different values 
_ of the total expressidn are found. The even values of 7 + 7’ and 
r — r’ produce as many more different values, hence the fractional 
differential coefficient of each of the een COS & e and sin x has 

Qn* different values.* 


For example, the differential iliac to the index + of cosz 
and sin 7 have 2 xX 2* or 8, values. Those of the former will be 
found to be : | 


Sr 
+ cos +2), +¥7— cos (= +2), 


+ cos (= + 2), cos (= $2). 


: The reader may. observe that if differentiation to the index + be 
performed twice in succession upon cos 2% ‘by any one of these 
formule, provided the same be used both times, the result will 
always be —sin 2, as it should be. : 


7. The same multiplicity of valuen attends the differential co- 
efficients of other trigonometrical functions. For instance, let us 
consider the functions e«”* cos nx and e”* sin mx. It will be con- 
venient to asume m = a cos 0, and =asin@. We may then 
change ax into z, and they become | 


e7 cos (a sin 8), °° % sin (2 sin 
which - we shall call © and 6’ respectively, as Mr. R. Murphy ee 
* M. Liouville did not employ the direct process which we have, and 


consequently did not discover that the differential coefficients had more than 
2n values. | 


| 
| 
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done in the Camb. Phil. Trans. vol. y, where he has proved some 
properties of them. We have then 

(x sin 0) + sin (x sin 0) 
(C08 9); therefore 
d* 
4 l = (con 0+ sin (009+ 
= $cos a (Qrm 6) +V—Ilsina + 
(x sin 0) + W— 1 sin (@ sin 0); 


cos $cos (2 sin 9 + Qarr + a) + sin sin 6 + 2arm + 
In like manner | 


= 


eT 005 8 (x sin -+ + a0) — sin (asin @ 4 Qar'm 4 
By addition and subtraction — 


d*Q 
Scos (r — r') ax + W— 1 sin (r — 7) and 
(V). 
| 


(r — aw + V— 1 sin (r — 1’) | 
08 sin + (7 + 7) ar + af}, 


As before, the number of values of each of these expressions 1s 
twice the square of the denominator of a, supposed a rational frac- 
tion in its lowest terms. | 


The necessity of supposing rand r’ different, may be proved as 
follows. If @ be changed into — 6, © remains the same and @' is 
altered only in sign: no new value of the differential coefficient 
‘ought therefore to be found by making this change: but that would 
be the case if 7 and 7’ had been supposed equal: therefore that 
supposition would not give all the values of the differential coeffi- 
cients. A similar remark applies to the last article. 


8. We now come y another division of our salaeek | in which we 
shall prove certain formule discovered by M. Liouville, and ap- 
plied by him to the solution of a variety of eee The first 


of these is 


foc + a) a™'da = (— (W). 


This may readily be proved by the theory of generating functions. 


> 
j 
Q 
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Let G (2) denote the generating function of g(x), or that 
function of ¢ in the expansion of which the coefficient of ¢ is 
(): then we have 


+a)= €*G g(a), | 
G/¢ (x + a) a""'da = fre a"-'da.G 9 (2). 
0 0 


Let e= £9, then a =i and the second member of the 


ens equation becomes 


(log 
but (log Go(2) =G (x): 


and + Ida = (— 1)" (2). 


It may be remarked, that if the form of @ (x) be such that the 
first side is infinite, the second may be made to agree with it by the 
aid of the complimentary function, if it do not without ; so that the 
formula is true for all forms of ¢(2). 


9. An obvious use of this formula is to find the form of ¢(z) 
_ when the value of the definite integral is known, as in the following 
problem. 


Each particle M of a line AB, infinite in both directions, exerts 
upon a particle P without it a force perpendicular to the plane 
_ABP, and proportional to sin PMA, and to an unknown function 
F (7) of the distance PM or 7. The total action of the line upon P 
is a known function f(y) of the perpendicular distance y of P from — 
_ AB. It is required to find the function F (r). 


This is an important physical problem, for if AB represent a 
voltaic current, and P a pole of a magnet, the action of M on P 
is such as has been described; and the total action is found by 
experiment to be inversely as y. 


Let C be the foot of the perpendicular from P on AB, aa let 
CM =z, then 7? = y*? + 2*, and sin PMA = ~ ‘whence the 
whole action of AB is | : 

y de FV 
| Vy? + Vy? +4 


which is equal to f(y) by hypothesis. 


* 


0 
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F 
Assume = 9(7*), therefore. 


+ 22) de 


Let 2, =a, then the first becomes 
+ a) am da,. 
which by is to 


Supposing that the differentiation is performed, and 7? is substi- 
tuted for xz, F(r), which is to is known. 


If f(y) = and by formula (E), p. 14, we find 
= 


a 


therefore F(r) = — 


It is easy to verify this result onan integration, and thus : 
to confirm the truth of the principles of the new calculus. — 


10. It will be convenient to possess a formula in which the 
limits of the integral shall be 0 and 1. Such a one may de ob- 
tained from Bal by changing the independent variable. 


Let etaaz; ; then when a=0, 0=1, and when a=o, 0=0: 


2% da= , and the is changed into 


also, a = 


| 
(1 — = (— 1)* (n) (2). 
Let aig (2) = (; therefore 


The following modification of the formula is often convenient. 


Assume = fa and = , then we have 


| 
> 
| 
4 


| 
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The differentiation indicated in the second member must be per- 


formed with respect to , as the independent variable. 


11. Asan example of the use of. the last formula, we select the 


following problem, which includes, as a particular case, that of find- 
ing the tantochronous curve. 


To find the curve such that the time occupied the of 


_a-particle sliding along it by the action of gravity to a given point 


from a given vertical height above it, shall be a given function of 
that height. 


Let 2 be measured vertically. Then the known expression for 
_ the time of descent —_ a height A to the origin, is 


ds 
on 


but by hypothesis this is equal to the known function ay. (h). 


Putting = = and the integral with that in 
formula (Y), we find | 


whence, observing that = 


By changing h into x, y (x) is known, and thus we mare | a diffe- 
rential equation to the curve. 


Suppose, as a particular case, that f (h) = ae then, after 


changing h into - - for convenience of one we have 


(n + 1) 


(n+) 
by (F), if 2 be positive. Therefore the equation to the curve is 
| | 
dx (n + 4) 


if n =, or the time be independent of the height, we have 


which belongs to a cycloid. 


| 
2 
| 
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If n = 4, or the time vary as the square root of the height, 


which belongs to a straight line. 
If m = 1, or the time vary as the height, 
2. 
dx 


putting ¢ ———" = = , and solving the equation, we obtain 


We might give many more instances of the use of these formule, 
but our limits will not allow us: and therefore we recommend to — 
such of our readers as are particularly interested in this subject, 
the original Memoirs above referred to, which are more complete 
in the applications of the Calculus than in its principles. 


I. 


VI.—ON THE TANGENT PLANE OF THE ELLIPSOID.* 


In an article, which appeared in the first Number of this Journal, 
some examples were given of the advantage of employing, in cer- 
tain cases, a form of the equation to the tangent of the ellipse, 
which does not involve the coordinates of the point of contact. 
There is an analogous form of the equation to the tangent plane to 
the ellipsoid, viz. | 
le + my + nz = + + nic’, 

1, m, n being the cosines of the angles which a perpendicular on 
the tangent plane makes with the axes. We shall not give a proof 
of this, but refer the reader to Hymers’ Geometry of Three Di- 
mensions, Art. 29. A few examples of its application may per- 
haps be useful. | | | 


1. To find the locus of the intersection of three tangent planes 
to an ellipsoid, which are mutually at right angles. - 


Let (lmn), (Um'n'), (l'm'n") be the cosines of the angles which 


* From a Correspondent. 
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perpendiculars on the planes make with the axes. Their equations 
will be | | 


4 my nz 4 mb +4 n2c?, 


| x 4 m'y 4 = pre re 
Squaring both sides of each of these equations, adding, and 
observing that, since the planes are at right angles, 
lm + + Um" =0, In4 +0 == mn + mn’ +m 
we have | 


This proposition was first proved by Monge. 


2. To find the locus of intersection of the perpendicular on the 
tangent plane with that plane. 


The equation to the tangent plane being | 
le + my + nz = 4+ mh? + (1), 
the equations to the perpendicular will be 


Multiplying the successive terms on both sales by these equal 
quantities respectively, we have | 


+ 22 = 4 + 
or (a? +y? 4 22)? = 4 by? 4 
3. Three planes, mutually at right angles, touch three eee 


whose radii are a, a’, a’, respectively. To find the locus of inter-— 
section of the planes. 


In this case the equations will be 
lx + my + nz 
lx + my + nz 
+ my + n'z =a". 
Squaring, and adding, we shall have _ 
4 2? — 4 a? 4 "2. 


4. Three planes, mutually at right angles, touch three concentric 
ellipsoids, whose principal sections have their foci coincident. To 
find the locus of intersection of these planes. 


Let abc, abc’, a’b'c", be the semi-axes of the three ellipsoids ; 
then the equations will be | 


|| 
§ 


4 
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lx + my + nz = + mb? + 
Ve + my + = + mp2 + n2¢2, 
+ my + = + + 
Squaring both sides, and adding, we have on the first side 
+ +2, as before. Also, 
[2q2 4+ m2b2 + — az 4 m2 (b? a?) n2 (c? 
19 + 2p'2 4 n'2¢c'2 — 4 [2 (a'2 b'2) 4 n'2 . 
: — since the foci in the principal sections coincide, | 
= 6% + (a? — b*) + n? — 
Similarly, | | 
4. mlb"? 4. — 4 (a? — + (b? — 
Therefore the second side becomes, by addition, 


= ai + + 
+ (I? 4. yor 
+ (m* 4. ml” n*) 


Similarly, | 
Hence the equation to the locus is 
gt + 22 =a? + 6% 4+ c” 


= q? — §%, 
or a? + 2 — a? + B, 
and b? + = 
therefore, by ies, 
And, in the same e way, | 
a + b2 4 
Consequently ‘ equation may be put into the fase 
xt + +b? 45% 4c") 
- This result was first published in the 19th volume of the 
Annales de Mathématiques. 


If the equation to the surface » 
= y = 


the caution to the tangent sei may we exhibited in the form 
1 (lx + my + nz) + 3 (am? + adn?) = 0. 
As there is no proof of this in Hymers’ Analytical Geometry, 
we subjoin one. 


“Again, since 


‘al 
j 
4 
a? 
‘ ‘ 
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_iet the equation to the tangent plane be 
le + my + nz = 6. 
This must be identical with the common equation to the tangent 
- plane drawn at a point (2'y’2'), viz. | 


Consequently the equation is 
(le + my + nz) + i (ma + na’) = 

The reader will find no difficulty in proving, by means of this 
equation, that tangent planes, mutually at right angles, will always 
intersect in a plane whose equation is | 

4 4 (a. + a) == 0; 
and that the locus of intersection of the tangent plane and the per- 
pendicular upon. it from the vertex, is a surface, represented by the 


equation 


+ y? + 22) +1 (ay’? + a2*) = 0. 


ax aX Xx 
: therefore we must have : 
| rz =-—G, 
| a! 
/ 
ax ao | 
| ma 
“whence 7’ = — —. 
Al | 
Similarly, 2? = — —. | 
| 
a a 
l 4[? 
2 
or § = — — (ma + na’). 
Al 
| 


12] 


_VIIL— TRANSFORMATION FROM RECTANGULAR TO. 
POLAR COORDINATES IN ‘DIFFERENTIAL EXPRES- 
SIONS. 


Wuewn we have a differential expression involving two independent 
variables, which we wish to transform into one with two other 
independent variables, the method to be pursued is not at. first sight 
obvious. If, for instance, we have the double integral | 


whites Vv; is a finetion of x and y, and we wish to transform it sO as 


to involve rand 0, the four quantities x, y, r, 0, being connected by 
the equations 
= f(r, = F(r, 8), 


dx dx dy dy 
ds = = dr +5, db, = + do, 
we cannot multiply the expressions for dx and dy jlo as they 
stand ; because in the integral y is supposed to be constant when x 
varies, and x to be constant when y varies. We must introduce 
one of these conditions by supposing dz = 0 when y varies, or 
=0 when varies. 


is the first of these, we have the equations 
| 


that 


O— 
dr da 

— dé. 


Eliminating da between these, we find 


do da dr dr dg 
_ From this it follows, that when dy = 0, dr =0. Hence 
= OD, 
do 
Substituting these values in the double integral, it becomes 


dxdy 
where V involves only 7 and 6. | 
When z and y are the rectangular, 7, @ the polar coordinates of 
a point, x = rsin 0, y = 7 cos 0, 
and da dy = rdr dé. 


If there be three independent variables the same method is to + he 
pursued ; but as we have to eliminate between three equations it 
becomes very complicated, even though we avail ourselves of the 
method of elimination by cross multiplication given in our first 


4 
| 
4 
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Number. In the particular case of the quantities representing the 
coordinates in space of a point, the assumption of a subsidiary 
quantity much facilitates the calculation. For let 


JIN dx dy dz 
be the function to be transformed, and 
cos 8, y =rsin@ sin sin 6 cosq. 
Assume p = 7 sin 0. Then we have | 
y=psing, =p sb), 


We can then transform, as we did before, from 2, y; z, to 
x, p, and after that from 2, p, @ to 7, 6, ¢. | 
The simplicity of the method consists in this, that the second — 
operation is exactly similar to the first, so that we do not require 


to repeat the calculation, but merely to write down the result, 


substituting r for p, and @ for Effecting the first operation, 


dy d 


dx dy dz = p da dp do. 
Again, changing from 2, py ‘te r, 0, we find 
dx dy dz = pr dr dp = dr sin 6 dé dg. 


dx dy dz = Vr dr sin dO dy. 
The same assumption may be usefully applied to transforming 
the expression 


Therefore 


av av 
dx? ar dz 


from rectangular to polar coordinates; for it will be seen, that 
9 


having found and Pr the whole expression may be written 


dy 
down without further trouble. 
Thus we find 
dy dp p dy 
dy? } dp* p? dg p dp 


4+ 


dp” dp) 


The expression for tied is got by putting 90 — @ f . h 
dz g ® tor @ in the 
above, so that without writing it down ee 
av 1 atv aV 


| 


dy 
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Putting rv for p and 0 for ¢, we have the similar expression 
dV dV dV 1dV dv 


| 
And the expression for -— in r and @ being similar to that for 


in 7 and 9, gives 
1dV_.1dV_. cot dV 
pe 

Adding these three expressions, 
da ap * (aw + 
1_@v adv 
r’sin?@ dg? dr* r dr 

r? d cos 0 Fant + GF’ 


which is the well known expression. 


IX._ON THE SOLUTION OF PARTIAL DIFFERENTIAL 
EQUATIONS. 


THE integration of Partial Differential Equations is much facili- 
tated by the principle which we have developed in our preceding 
Numbers, and it is the more remarkable, that it has been so little 
applied to these equations, as the first step was taken many years 
ago.. Fourier, in his 7raité de Chaleur, published in 1822, has 
shown that the series which are obtained in the solution of several 
partial differential equations may be conveniently expressed by the 
separation of the symbols of operation from those of quantity. But 
though he has used this method very frequently, yet he appears to 
have had some unwillingness to give himself up to it entirely as a 
guide in his investigations, as if he were not familiar with the prin- 
ciples on which it is founded. His idea apparently was, that the 
expression which he obtained as solutions might be conveniently 
expressed by separating the symbols of operation, and not that the 
symbolical expressions are the proper solutions of the equations, 
and the series merely the expansion of them. Other French 
writers seem to have avoided carefully entering at all on the track 


| 
2, 
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which Fourier opened: Poisson in particular, in a digression on the ~ 
subject of partial differential equations in the second volume of 

his Mécanique, does not put in the symbolical form the solution of | 
a very simple equation, which is so given by Fourier. Mr. Great- 

_ heed, in a paper published in the number of the Philosophical 
Magazine for September 1837, was, we believe, the first to call the 
attention of mathematicians to the utility of this method in the case 
of partial differential equations, but he had not then reduced it to 
its greatest degree of simplicity ; and his paper is chiefly occupied 
with a particular class of equations of the first order with variable 
coefficients, which are not so interesting as many others with con- 
stant coefficients. We shall here, therefore, proceed to give several 


_ examples of the application of the principles which we laid down 
in Art. V. of our first Number. 


And, first, we may observe generally, that linear partial differen 
tial equations between any number of variables with constant coef- 
- ficients, are to be treated exactly like ordinary differential equations 
with regard to one of the variables, the symbols of operation of 


the others being treated as constants. If, for instance, we have 
the equation | | 


This may be put under the form 
d 
(a +65) z= 0, 
| d | 
+05) (c + 0), 


where we suppose 2 to be the variable, and 7, * constant with 


and therefore 


regard to it. Now the operation (« = +b =) by the theorem 
given in page 25, is equivalent to | | | 
| | | 
a dy [dx dy | 
which being performed, gives 


S (y) being an arbitrary function of y, taking the place of the 
constant in ordinary differential equations. 


From this expression we Bet 


= — bx), 


dz dz 
~ 
+ b6—=c 
d 
~ 
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as by Taylor's theorem 


h. — 
e f(x) = f(x +h). 
We might equally well have supposed y to be the variable, and £ 
a constant with regard to it. 


Again, taking the well-known equation for then motion of waves, 
d*z d*z | 


d# 
it may be put under the form 


d? —d2\ 


and integrating it like the ordinary differential equation 
| | 72 | 
— 
| at —at 
we find z=« 


or + at) + — at). | 
The equation r — 2a s+ a*t=0 may be put under the form _ 


the solution of which } is by the Theorem in page 25, 


W f2 du2.0 


and functions of y arising from the 1 in- 

_ tegration. Hence, finally, we have 


z= ap(yt ly + ax). 
The equation 


r — at 4 2abp + = 0 


is equivalent to 


( — sal)! (2 


Integrating with regard to the first me we have 
d d ala 


d d ) 0 | 
— —a—)z=0, 
dx d 
Integrating with regard to the remaining factor 
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d 


| 


+ — ax) 
by changing the arbitrary function _ 
— tax) (y + ax) + — ax) 
= + ar) + — ax). 
Let us take also the equation 
| LY. 
Without operating on each factor separately, we may arrive more 


readily at the result by the same means as those employed in 
page 28 of our first Number. For we have : 


= Cv) + oy + — ax) 
¢ dy (xy) + o(y + ax) — az) 
d2? d-2\—' say | 


Therefore 
+ + — az), 

the arbitrary functions being added as in the second example. 

A very simple equation being one which occurs in the theory of _ 


heat is 
dy 
dt 


which is the expression for the rectilinear propagation of heat. 
The solution of this is easily seen to be, if we integrate with regard 
| | 
T= das I(x), 
f(z) being an arbitrary function. If the sign of operation be 
expanded, we shall obtain a series which is the solution derived by 
Poisson from the method of indeterminate coefficients. Laplace 
has deduced from the series an elegant expression for the solution 


under the form of a definite integral; but it may be more easily 
deduced from the symbolical solution in the following manner. 


Since = Vz. 
and also b? — WV 
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we can put the expression for v under the form 


= fix + Vat) 


by Taylor's theorem ; and this is Laplace's expression. 


2» 
This equation o = =a = not being homogencous i in the index of 


the operations, admits of two solutions of very different characters. ? 
The one, which we have found by integrating with regard to ¢, con- 
tains only one arbitrary function of x. The other, which may be 
found by integration with regard to 2, must contain two arbitrary 
functions of ¢, as the index of operation is of the second degree. If 
we write the equation in the form | 


dv ] 


and integrate by the method employed in page 28, we find for _ 
the integral 
(1 +2 ia 1.2.3.4 &e.) 


It seems at first sight anomalous, that the same equation should 
have two solutions so different in character: the following is the 
explanation of the difficulty. Since by Maclaurin’s theorem any 
function of a variable may be expressed by means of its differential 
coefficients taken with regard to that variable, for the particular 
~ value O of the variable, we know the function if we can determine 
its successive differential coethicients. Now from the equation © 


dv es oe we can, when we know the value of v when t=0, 


dt dx? 

determine the values of all the differential coefficients with 
regard to ¢ when ¢=0. So that in the resulting expression 
deduced from Maclaurin’s theorem there is only one quantity left 


undetermined, which is the arbitrary function f(x), introduced in the 


d% 1 dv 
integration. But from the equation a? We can only 


determine from the value of when x = O the values of the alter- 
nate differential coefficients with regard to 2. There must conse- 
quently be introduced another undetermined quantity, namely, the 


v | 
value of = when x = 0, for knowing these two quantities we can 
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determine the values of all the successive differential coefficients 
with regard to 2. : 


equation for the motion an elastic 
3 


_ The solution of which is readily seen to be | 


v = cos os (¢ F(z) + sin 


The equation for determining the vibratory motion of elastic — 
plates is | 
d%y diy d4y 
dat dy? dy* = 0, 
which may be put under the form | 
dv d? 
the integral of which is 


v= cot (7 + +S) + 


In investigating the motion of heat in a ring, we gente as an 
equation for determining the temperature at any time and point, 
dv d*y 
dt 


the solution of which is 


— 


an expression ‘dead connected with one previously given. 
In the examples which we have given, the coefficients are con- 
_ stant, but if one of the variables only enters into the coefficients, 
and we integrate with regard to the other, as the one variable is 


unaffected by the sign of operation with regard to the other, it may 
be considered as a constant in the integration. 


A good example of this kind of equation is that which expresses 
the motion of heat in a solid cylinder of infinite length, namely 


dv 2a I] =) 
the solution of which is — 


at +2 f(2), 


by 
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or, as it may be expressed, 


But this will only be possible when the integration is of the 


first degree, since, as we showed in our first Number, the sym- _ 


bols of operation become subject to different laws when the va- 
riable itself and the sign of differentiation are both involved. 
And this leads us to the consideration of equations with variable — 
coefficients. As in the case of ordinary differential equations, the 
solution of this class is attended with great difficulties, so as to 
become almost impossible for equations of an order higher than ~ 
_ the first, it is to be supposed that these difficulties are no way _ 
diminished in the case of partial differential equations. Of these, — 
however, when of the first order, Mr. Greatheed has shown, that 

a large class may be solved like ordinary differential equations. — 

It is included under Py form 


+e 


where X is a function of x . Y of y only, and P and Q 
of both x and y. gg _ be reduced to the form | 


ar +a 
dy dz _ dz 
then Y a7 


so that it is only necessary to consider’ the latter equation. | Let 
it be put under the form. 


ly 


and treated as an ordinary linear differential equation between 
z and x. If we integrate by the method of integrating factors, 


by a change of the variable; for if y' = 


| d 
the factor is (x and the 


whence | | 
dr(X — — 
( (y) being an arbitrary of y. 
We shall take as an ae the equation 


= dt, and therefore y= et, “Then 
de ' 2 


| 
* 
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The ti of which, by the — formula, is 

or z gn log ¢ (t) 


log x) = 2" y ~ log2) = any (2). 


may be integrated in the same way. Or if we ‘diange both the 
dy 


independent variables, making a dt and a = du, it becomes 
(= x) 
which gives z= ¢ “/du dt + e “4 (t) 


d 


zy 
as u = log az, t = log y. 
The equation x 
qu | + dy 


a8 may, by changing sg into 4 d.y’, be transformed into 


d ] 
+ [2 ay z= 
The of is | 


2—€ 


I 


(x + y) p(y? — 2°). 


| 


3 
¢ 
—7? —— 
o(y*) 
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Other equations, which at first sight do not appear to come under 
this form, may be reduced to it by a proper assumption of a new 
independent variable. instance, the 


is converted into | | 
dz 9 dz 
by supposing y =u — 2. These however are particular cases 
which come under no > general rule. | 


F. G. 


X.—ON STABLE AND UNSTABLE EQUILIBRIUM.* 


Ler P, P’, P”...... be any sauiber of forces acting upon a material 
system, and let dp, dp’, dp’...... be the virtual velocities of their — 
points of application, estimated positive in the directions in which 
the forces act; then, when P, P’, P’...... are in equilibrium, the 
equilibrium is stable when S(Pdp + Pdp' + ...... .) is a maximum ; 
and unstable when this is not a maximum: and vice versd. | 

We shall demonstrate this in the case of one rigid body; and the — 
same reasoning can easily be applied when the system contains 
several rigid bodies. | 

When P, P’, P”...... act upon a rigid body, they can always i 
reduced to two forces, but not to one. Let R and R’ be these 
forces, and dr, dr’ the virtual velocities of their points of applica- 
tion: then | 
Pdp + P’dp’ R'dr’, 
= {(Pdp 4+ P’dp’ ) = {(Rdr + R'dr’). 


Hence du — Rdr + R'dr, 
and + dRdr + R'd?r’ + dR’dr’. 
Now when P, P’, P’...... are in equilibrium, 
Pag + + ...... = 6 


Rdr + =0, ..du— 0; 
and, since dr ae dr are independent of each other, 
R = 0 and R' = QO, 
du = dR.dr 4 dR'.dr’. 
When wu is a maximum, du is negative, and therefore (since dr 
and dr' are independent of each other) dR.dr and dR’.dr’ are both 


* From a Correspondent. 


| 
Z 
f 
f : 
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true. 
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negative. Hence if dr be positive, dR is negative; and vice versa: 
and so with dr’ and dR’. 


_ From this we learn, that if the body receive any slight displace- 
ment from the situation of equilibrium, two small forces dR and dR’ 
are brought into play, which act opposite to the directions in which 
the points, at which they act, have moved in consequence of the dis- 
_ placement of the body. Hence the effect of these forces is to draw 
the body back towards its position of oe and meee the 
equilibrium was stable. 


Conversely, i in order that the equilibrium may be stable, the forces 
dR, dR’ put in play by the displacement must act so as to carry the 
body back to its situation of equilibrium : and therefore the virtual 
velocities dr and.dr’ must have different signs to dR and dR’ respec- 
tively: hence d*u is negative, and therefore u a maximum. 


If w be a minimum we may shew in the same manner that the 
equ ilibrium is unstable, and vice versa. 


If u be neither a maximum nor a minimum, then d?u (or d°u, if 
du vanishes) will be positive for some displacements, and negative 
for others; and will therefore be apparently dubious, though in fact 
it will be unstable ; for although a displacement may be found which 
shall cause the body to return to its position of rest, yet in oscillat- 
ing through the position of rest it may come to one, from which it 
will not return to its position of rest. The converse of this is also 


Hence, if by stable equilibrium we mean, that for all small dis- 


_ placements the body oscillates about its position of rest, then the 


equilibrium is stable when w is a maximum, and in no other case: 
in all other cases the equilibrium is unstable. 

N.B. ~~ We have supposed that Pdp + P’dp’+... isa perfect dif- 
ferential. If this be not the case, the above reasoning is still true ; 
but the result should be stated thus: when the equilibrium is stable, — 
the first differential of Pdp + P'dp'+... (or if that and the second 
- yanish, the third) i is negative, and vice versa: in all other cases the 
equilibrium is unstable. | 

If the system consist of several bodies, and R, R’, R’, R”’...... 
be the smallest number of forces to which the system of forces act- 
ing on the bodies can be reduced, then 


du = Rdr + R'dr + R'dr" + + .. 


in which dr, dr’, dr’, dr’” ...... are independent of each other: and 
the reasoning upon these and the forces R, R’, R", R" ace will be 
as above. 


The following is an application of this ‘oiloaiule to fluid equi- — 
librium. 


PROBLEM. Dinies the form of equilibrium of a mass of | 
incompressible fluid, every particle of which is attracted towards 
a centre of force, varying as a function of the distance; and to 
determine whether the equilibrium be stable or unstable. 
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Let the fluid be referred to polar coordinates, the centre of 
force being the origin, @ the angle which r makes with a fixed line 
through the origin, and w the angle which the plane through r 
and the fixed line makes with a fixed plane passing through the 
fixed line: and let (as usual) cos @ =p. It is evident that the 
fluid must completely surround the centre of force: let r andr’ 
be the radii of the external and internal surfaces; these are 
functions of » and w. The volume of the fluid is given: let 
c be the radius of the sphere which has the same volume: hence 


SF dw dr = 9) du dw. 


_ Since r° and r? are functions of p and w, they may be expanded in in 
series of Laplace’s coefficients: let 


r= ai4a(Y, + Y, + ay, 
and — a3 4 ) = a? + ay, 


_ where a is a numerical coefficient, the use of which will wee seen 
hereafter. 


Then, by a fundamental property of shine functions, 


SL y tude =0; 
and similarly of 7/. po we have, by the above equation, 
i S(r) be the force of attraction on a unit of mass at a 
distance r from the centre of force: then f(r) em is the attraction 


on the element of the mass ém:: this attraction acts in the line 7, and | 
tends to shorten 7: hence, in this case, 


u= [%.3— f(r) dm? dr; and this = 3./$— f(r) dr} dm. 

The symbol of integration = refers to the differential of the 

mass 6m. Now ém equals (as above, if density = 1) 
sin dw dr = — r* du dw dr, 


and the limits of 7, w, p are respectively r and r, O and Qn, 
1 and —1: then, putting 


f—f(r) dr =fi (7); and dr = 


we have 


72 f(r) dy dw dr = $9 (13) — du dw. 


But 9 (1°) — 9 = g(a + ay) — + ay) 

= 9(a°) — + ay (a*) — ay (a?) 
(0°) — (a) +... 
in which the accents to @ denote differentiation with respect to 
a’ and Hence 


= (a?) — 40 (a3) 
+ (a2) — dy de +. 


> 
| 


“184 Analytical Geometry of Three Dimensions. 


It is very easily shewn, that (a°) = a5 (a), and is there- 


fore negative in every case of attraction, (but positive if the force 
were repulsive ). 


From these calculations we gather the following results. If a be 
a very small quantity, the increment of wu above what it becomes 
when a = 0 (that is, when r =a and r' = a’) does not involve 
the first power of a: hence, when r =a and r’ =a, we have : 
du = 0, and therefore the fluid is in equilibrium. — 


Hence a form of equilibrium is that of a hollow sphere of any 
dimensions (so long as the volume be constant), and having the 
centres of the bounding surfaces at the centre of force. 


And the equilibrium of the ezternal surface is stable, because 

for all values of y when y'=0 the multiplier of a? is negative, and 
therefore uw is a maximum: and the equilibrium of the internal 
_ surface is unstable, because for all values of y’ when y = 0 the 
multiplier of a’ is positive, and therefore u is a minimum. On | 
the whole, the equilibrium i is dubious. | 


Cor. 1. If a'=0, then /™ = 0, and the equilibrium i is ge 
stable. 


| Cor. 2. If the force were pen rae fee the fluid were held 
“sae, in equilibrium hy a rigid spherical envelope; then y=9%, and the | 
| multiplier of a* would be negative, since g” (a*) is then positive, — 
and u is a maximum, and the equilibrium of the internal surface i is 


stable. 


P. 


XI. —ANALYTICAL GEOMETRY OF THREE DIMENSIONS. 
No. II. 


WE proceed to give some instances of the advantage of symmetry 
in applying the Differential Calculus to Geometry of three dimen- 
sions. The following form of the equation to the tangent plane is 
probably known to many of our readers; but, as Leroy has only 
deduced it from that in terms of the partial differential coefficients — 
of one of the coordinates with respect to the other two, we will give 
an independent and easy proof of it. | 


1. If F(a, y, z) = 0, be the equation to a surface, the locus of 
the tangent lines drawn to it at a point (2, y, z) is a plane whose 
equation is 


© 
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equation to any line through the point x, Zz, is 
| m n 


other points. Let the coordinates of the nearest of these be yp Yp 2 
and let each member of equation (2), when 2’, y’, z’ become 2; Yy 2p 
be assumed equal to 7, then 


y=y + mr, 
F(a), = 9; 
therefore F lr, y+mr, z4+nr = 0. 
Expanding and observing that F(x, y,z) = 0, 


When the a tangent, the pointe (x,y; 2)> (Xp Yv Zp) 
approach indefinitely near to one another, and 7 becomes indefinitely 
small ; therefore, taking the limit of the preceding equation, 

dy 

Multiplying the several terms of this by the several members of 

equation (2), we obtain 


+ &e.) 4 Oe 


for the equation to the locus of the tangent lines. 


2. Since the normal line is perpendicular to the tangent plane, 
its equations | 


dx dy dz | 
3. We shall next investigate an expression in terms of the partial - 
differential coefficients of F (a, y, z), for the radius of curvature of 


a section of a surface made by any plane passing through the nor- 
mal at a given point. 


Let = U, dy = V, W, then the equations 


‘normal are 
Again, if ds be an element of the section, and dz, dy, dz its 
2 
ds’ ds’ ds 


projections, are the cosines of the angles which it 


- This line will in general be cut by the surface in one or more ~— 


; 


| 
/ / 
; 
* 
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makes with the axes, and therefore the equation to the normal | 
plane to the section is 


(a! 2) de + y) dy 4 — 2) de =0 (8). 


The line of intersection of two consecutive normal planes will be — 
determined by (5) and its differential, which is ) 


+ (y —y) Py + —2) —ds* = 0......(6); 
and the intersection of this line with the normal to the surface will 

evidently be the centre of curvature of the section. We may sup- 

pose 2’, y’, 2 to belong to this point; they will then be the same in 

(4) as in (5) and (6): and assuming each member of (4) _ to 

Q, and substituting for x’, y', z' in (6), we have 

(Oda + + Wadz) Q — ds? = 0. 

_ But if p be the radius of curvature 


therefore 


+ V2 4 W?.ds? 

= = + Vdy 4 Wd 

This expression may be transformed as follows. Since 
U dx + Vdy + Wdz = 0, 


and dy? =% =, 
dzdx” dedy | 

U 4+ Vdty + 4+ + + wdz 
+ 2u dy dz + Qv' dz dx + 2w dx dy = 0. 
Employing this equation, and assuming = I, 2 = m, —n, 
(7) is changed into ; | | 
| V7U2 + V2 4+ W2 


Py + my + n2w + 4+ + 2lmw' 


4. It may be as well here to show how formule involving the 
partial differential coefficients of F, may readily be transformed into 
others in terms of the partial differential coefficients of one of the 
coordinates as z. Suppose the equation to be put under the form 
0, so that F (a2, y, z) = f(a, y) —z. Then 


d¥ dz dF dz 
dr dy dy’ dz 
dt d% @F 
(9), 


dydz dzdx dady  dxdy’) 


— 
7 
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If we substitute these values in (8) and adopt the usual notation, 
we obtain the known expression 


Vit 
p= > 
Pr + + mee 


5. To find the ‘tint whole radius of curvature is a maximum 
or minimum, and the values of those radii. : 

In formula (8) let + Vv U* + V* + W? =P, then 
Pp 


the quantities J, m, n being connected by the two equations 
+m+nt =1. 
1U + mV +2W =0 
Differentiating these equations, 
(lu 4+ mw’ + nv’) dl + (lw + mv + nu)dm 


+ (lv ml + nw) dn = 0 
Udl + Vdm + Wdn = 0.... .... (15), 


(13) + (14) x A+ (15) gives, on to zero, the coeffi- 
cients of each differential, | 


lu + mw + nv + NM + pU = 0 
lw’ + mv + nu’ + Xm + = (16) 
lv’ + mu + nw -+ AN + pW - = 0 


Multiplying the first of these by J, the second by m, and the 
third ay n, adding, and — by (10), (11), (12), 


A == 0. 
Substituting for \ in (16), they become - 
(uw —=) 04 wm -+- vn + = 0 
| p | | 
/ 
p | 
vl + um (w — n 


To obtain the equation for p, it remains to diisideke l,m, n, yt 
from whe coe (17) and (12). The result is 


tq 
. 
ef 
< 
. 
é 
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— ( --) — 2WUe' (« 2UVe ( 


p 
2VWo'w’ + 2WUw'w' + 2UVu'e' 


6. If the equation to the surface ba of dhe rue 
o(z) +x(y) + =0, 


which includes, among others, all the surfaces of the second order 
when referred to their principal diametral planes ; U, V, W. are 
functions of x, y, z alone, respectively, so that w’, v', w’ are all 
zero, and the equation (18) simplifies considerably. In short, in 
this case, equations (17) give immediately : 


“—— w— — 
p p 
and by substituting these values in (12), we have 
U? W? 
p p | 
«7. ExAMPLE. To find the principal radii of curvature of an 
ellipsoid. | 
In this case ae 
2 2 
-1=0, 
On Qz 2 2 
and the value of P is a | 
xr 22 


if p ‘denote the perpendicular from the centre on the tangent plane. : 
Substituting these values in (19), we have. 


y* 
a (pp — a) (pp — * (pp — 
which may easily be put in the form | 


P 
p p 
| 
| 
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8. As another example, we propose to find the principal radii 
of curvature at any point of the surface whose equation is 


+ + = (2° + + ...... (a), 
and which is the locus of the extremity of the perpendicular from 
the centre on the tangent plane to an ellipsoid. Put 
a2 4 22 = 7, 
then F (a, y, = 4 by? + — rt, 

U 22 (a? — 2r*),, V =2y(b?— 2r?), W = 22(c? — | 

2 (a?—2r?— 427), v= 2 — 2r?—4y?), w= 2 (c?—2r?2— 422), 
w= — Byz, = — 822, w = — Bry. | 

Instead of substituting these values in equation (18), we shall 
employ the equations from which that is deduced, namely (12) and 
(17), because by this method it is more easy to obtain a result of 
-asimple form. Of these, (12) becomes 


lx (a? — Qr?) + my — + nz (ce? — 0 (8), 
and the first of (17 yy? 


2) — aay + px (a 298) = 0, 


or 1 G — 27? — ~) — 4x (la + my + nz) + pa (a? — 27?) =0. 


then we have | 


l G — Qr2 — — 4he + pa (a? — 2r?) 0 | 
m — 278 + py (8 — 272) = Ob 


— 2r? — -) — 4kz + pz 2r*) = 0 


and we shall obtain the equation for p by eliminating J, m, n, k, p, 
from the five equations (6), (c), (d). For this purpose, multiply 
equations (d) by 2, y, z respectively, and add them, reducing by 
(a), (6), then we find | 


+ 472) + pr’? = 0. 

By this relation, eliminate » from equations , and we have — 
4; (a2 — r*) 72+ (a? — 

a? — — 2 

with corresponding expressions for m and n. Substituting them in 
(6), and denoting 


F (a, x) + F (6, + F(o z) by 2). 


‘= 


| 

} ug 

4 

4 

« 
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4 (a? — 7?) (a 2r?) r2 (a? — Qr*)? 


| 
but by (a? — 1°) a? = 0: 
multiplying the last equation by 472, —— from the preceding, 


and dividing the result by — ait 
4 
— 


p 
| which i is the simplest form in which the equation ¢ can be presented. 
The value of P will be found to be | 


2 4 bay? e422, 


The of the perpendicular from the centre on the tangent 


plane to this surface is 


4 
therefore if p denote that line, P = =. 


9. The directions of the sections of greatest and least curvature 
will be found by eliminating \ and » from —— (16). The 
result is 


(lu 4 mw! + nv') (mW — nV) + (lw + mv 4 nw’) (nU — IW) 
+ (lr + mu’ 4- nw) (lV — mU) = 0, 

(Vol — Ww') 2 + (Ww' — Uw’) m? + — Vo’) 

—Vw' + U (v—w)} mn + — Wu' + V (w—u)} nl 
+ 3Vu'—Uo' 4 W (u—v)3 lm = (20), 


_ which equation, together with (11) and (12), determines the values 


of 4, m,n. It is easy to see that there will be two sets of values of 
m?, n*, and no more. 


10. To prove that the directions of greatest and least curvature 
are at right angles. 
For abbreviation, write the equation (20) 
+ Mm? 4+ Nx? + Limn + M'nl + Nim = 
Substituting the value of m derived from equation (12), 
(LE + Mm? + N'lm) W? — (Lim + M7) (WU 4+ mV) W 
+ N(lU 4+ mV)? = 
or (MW? + NV? — L'VW) m? + SN W (L'U + M’Vy ln 
+(NU24 LW?_MWU)2—0, 
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| | 
a quadratic equation in 73 whence we find, that if /,, /, are the 


two values of J, and m,, m, those of m, 
mm, _ NU? 4+ LW? — M’WU_ 
Ul, MW? + NV? — 
so that if we assume 
il, = K (MW? + NV’ —L'VW), 
shall have = K (NU* + LW’ — 
and also = K (LV? + MU? - N’ UV). | 
Hence, the of the angle between the two directions, 
Ll, 4+ + nn, = K $L + W?) + M(W? + U?) 
+ 'N (U2 + V?) — L'VW — M’WU — NUV}. 
The value of the quantity within the brackets is zero, as may 
easily be seen by recurring to equation (20) for the values of 


LM, &e, and collecting the terms multiplied by w, v’, &c. respec- 
--Therefore the ‘two directions are at right angles. 


To find the lines of curvature of a given surface, or the 
; Prior the in which we must pass from one point to a consecutive, in 
order that the normals at the two points may meet. 


Retaining the same notation, the equations to the normed at 
point x, are | 
| 
Let each member be assumed equal to Q, then 
a—=2x+QU, ¥=y7+ QV, 2 =2z + QW. 
If 2’, y’, 2’ belong to the point in which two consecutive normals 


meet, the differentials of these three equations must be verified at 
the same time with them, therefore | 


dy + QdV = 0, 


dz + QdW + WdQ = 0, 


or udx + widy + vide + de + 
wdx + vdy + 4. a+ + V = 0, 
vdx + + + wa 


These equations are of the same eae as equations (16), therefore 
the elimination of Q and dQ from these must give the same relation 
between dz, dy, dz, as the elimination of \ and p from (16) gave 
between /, m, n, namely (20). And if 


+ dy*® + dz*? = ds’, 
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=. 2, fi will be the cosines of the angles which an element 
of a line of curvature makes with the axes, and they satisfy the 
same system of equations as 4, m, n, namely (11), (12), (20). 
Therefore the directions in which two consecutive normals meet 
coincide with those of greatest and least curvature. 
For the sake of shortness, we shall eliminate Q and dQ from 
equations (21), instead of their expanded forms. This is easily 
done by cross-multiplication, and the result is 


U (dVdz — dWdy) + V (dWdx —dUdz) | 
+ W (dUdy — dVdx) = 0......... (22): 


from this equation and that to the surface, the differential to the 
projection of a line of curvature on ei coordinate plane may be 
found. 


12. Exampxe. To find the lines of curvature of an. ellipsoid. 
Let the equation be 


2x Qz 
then U=—, 


aU = av dy dW = 


Substituting i in (22), 
(0? — c?) x dydz + (c? — a*)y dz dx + (a? — zdxdy =0. 


Multiply by = and substitute the values 


ge y dy 
therefore 
dy\ 
| 2 
b*) (1 dx dy = 0, 
or, reducing, 
| 
ay 


| 
— — (a — de dy ry = 0. 


To integrate this, multiply by 4ay, and let x? = w, y* =», then 


(a* — du dv — 


u dv? 


(a? c? 


= 
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bt —c* a | 
or ( dv + du) (udv—v du) —(a’ — 6°) dudv=0 
du 
| 
au ) dv 
az — c2 du’ 
b? 
an 1 equation of Clairaut’s form, whose solution is 
| C (a? — 6?) 
| C (a? — 32) 


the equation to an ellipse or hyperbola. For the method of deter- 
mining the constant, see Leroy, No. 422. : | 


XIl._—_MATHEMATICAL NOTES. 


1. Taylor's Theorem. The following proof of Taylor's Theorem 
_ by means of the separation of symbols of operation from those of 
quantity, may be interesting to some of our readers. Let there be 
an operation, which we may for convenience call D, which per- 
formed h times on f(x) converts it into f(a +h); that is, let 
Di f(z) =f (x + we wish to determine the nature of D. 


definitely 


— f(z) D'—1 


when h is indefinitely small. clan taking the symbols of 
operation separate, : 


dx h 


Now we know that — when h is in- 


when h = QO, 


d\h 
which gives D = (1 th =). when h = 0. 


Expanding by the oe Theorem, and making h = 0, we find 


Ps 

‘id 

* 

rE 

if 

ay: | 
~ o- a 
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or D 


therefore f(x + h) = Sf (2), 


which is Taylor's Theorem. 


2. Problem from the Papers of 1835. If p, r be the perpen- 
dicular from the origin on the tangent plane and the radius vector — 


2 
of any surface, then . will be the perpendicular on the tangent 


plane at the corresponding point of the eater, which is the locus 
of the extremity of p. 


Let x, y, z be the coordinates of the first surface, a, P, y of the 
second. As p is the perpendicular on the tangent plane, the equa- 
tion to the plane is 


ax + Py + yz = p? = a? thay oseeee (1). 

Since this plane is a tangent to the surface, this equation will 
hold good if for x, y, z we put «+ dz, y + dy, z + dz. Whence 
adx + + ydz=0...... (2). 

Now if V =0 be the equation to the surface of which a, Ae 
are the coordinates, and P be the = on the tangent 


plane ; then 


| dV 


~ dv Vv 
aa) * (aa) + 
a¥ 
Now differentiating (1), ty [2 yz as 
and paying regard to (2), we have 


— 2a) da + (y — 28) dB + (z — 2y) (4), 
A (3)—-(4) gives, on equating to 0, the coefficients of each of the 
differentials 


ay... dV dV 


Substituting these in the expression for P, it becomes 
P = 2 (a? + + 7) — (on + 


which by (1) is reduced to 


+m 


\ 
‘ 
— 
da dij dy 
7i 
* 


